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ABSTRACT 

The calcium and hydrogen abundances, Galactic positions and kinematics of 146 DZ 
white dwarfs from the Sloan Digital Sky Survey are analyzed to constrain the possi- 
ble origin of their externally polluted atmospheres. There are no correlations found 
between their accreted calcium abundances and spatial - kinematical distributions 
relative to interstellar material. Furthermore, two thirds of the stars are currently lo- 
cated above the Galactic gas and dust layer, and their kinematics indicate multi-Myr 
residences in this region where interstellar material is virtually absent. 

Where detected, the hydrogen abundances for 37 DZA stars show little or no cor- 
relation with accreted calcium or spatial - kinematical distributions, though there is a 
general trend with cooling age. It is found that Eddington type accretion of interstel- 
lar hydrogen can reproduce the observed hydrogen abundances, yet simultaneously 
fails to account for calcium. The calcium-to-hydrogen ratios for the DZA stars are 
dominated by super-solar values, as are the lower limits for the remaining 109 DZ 
stars. All together, these polluted white dwarfs currently contain 10 20±2 g of calcium 
in their convective envelopes, commensurate with the masses of calcium inferred for 
large asteroids. 

A census of current T c g < 12 000 K, helium-rich stars from the Sloan Digital Sky 
Survey suggests the DZ and DC white dwarfs belong to the same stellar population, 
with similar basic atmospheric compositions, effective temperatures, spatial distribu- 
tions, and Galactic space velocities. Based on this result, pollution by the interstel- 
lar medium cannot simultaneously account for both the polluted and non-polluted 
sub-populations. Rather, it is probable that these white dwarfs are contaminated by 
circumstellar matter; the rocky remains of terrestrial planetary systems. 

In this picture, two predictions emerge. First, at least 3.5% of all white dwarfs 
harbor the remnants of terrestrial planetary systems; this is a concrete lower limit 
and the true fraction is almost certainly, and perhaps significantly, higher. Therefore, 
one can infer that at least 3.5% of main-sequence A- and F-type stars build terrestrial 
planets. Second, the DZA stars are externally polluted by both metals and hydrogen, 
and hence constrain the frequency and mass of water-rich, extrasolar planetesimals. 

Key words: circumstellar matter — minor planets, asteroids — stars: abundances — 
stars: chemically peculiar — stars: evolution — planetary systems — white dwarfs 



1 INTRODUCTION 



The curious presence of heavy elements in the atmospheres 
of cool white dwarfs has been known since the discovery of 
the first degener ate stars; van Maa nen 2 is the prototype 
DZ white dwarf llvan Maanen|[l917l ). For more than half a 
century it has been understood that, in the absence of ra- 
diative forces, elements heavier than hydrogen and helium 
should rapidly sink below the phot osphere in the hig h grav- 
ity environment of white dwarfs IjSchatzma ni l 19481 ): hence 



metal-lined white dwarfs like vMa 2 must be externally pol- 
luted. 

Owing to a relatively low atmospheric opacity, the de- 
tection of calcium in the optical spectra of more than one 
dozen cool, single, helium-rich (DZ) white dwarfs preceded 
the availability of 8 m class telescopes and high resolution 
CCD spectrographs (|Sion et al.lll99 0b). Although the sink- 
ing timescales are always orders of magnitude shorter than 
the white dwarf cooling age, heavy elements can persist for 
up to 10 6 yr in the sizable convec tion zones of helium-rich 
degenerates l|Paauette et al.lll986t) . This fact led somewhat 
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naturally to the hypothesis that the source of the accreted 
heavy elements was the interstellar medium (ISM), encoun- 
tered in sufficient density a few to several time s per Galactic 
orbit ( o n avera ge) for the Solar neighborhood l|Dupuis et al.l 
ll993al lbl [T993 ). 

Several factors have continually challenged the ISM pol- 
lution scenario. First, the typical dearth of hydrogen rel- 
ative to calcium in DZ st ars argues for the accretion of 
volat i le depleted mater ial l|Dufour et al.1 120071 ; IWolff et all 
120021 ; ISion et alj|l990ah . Second, the firm establishment of 
the hydrogen-rich DAZ spectral class, with heavy element 
diffusion timescales as short as a few days, and Galactic 
positions far from known inte rstellar clouds l|Koester et al.l 
l2005al ; IZuckerman et ai]|2003j ). Third, the increasing num- 
ber of metal-contaminated white dwarfs with infrared ex- 
cesses and closely o rbiting rings of dusty and gaseous debri s 
IIFarihi et al. I l2009l ; lJura et al.ll2009al ; iGansicke et~aIll200Sl ). 
Thus, while there is firm observational evidence of pollution 
via circumstellar material, as yet there is none favoring the 
ISM. 



lAannestad et all (| 1993T ) provided an extensive spatial, 
kinematical, and calcium abundance analysis against which 
they tested likely scenarios of interstellar accretion based on 
the best available data at that time, concluding interstellar 
accretion could not explain the abundances in any of their 
15 DZ white dwarfs. However, their landmark study was lim- 
ited to targets within 50 pc for the most part, and to stars 
that were proper-motio n selected, and he n ce biased toward 
higher space velocities l|Sion et al.l Il988h . iKilic fc Redfieldl 
120071 ) examined possible correlations between the current 
positions of 35 DAZ white dwarfs and models of the sur- 
rounding ISM as inferred from nearby column density mea- 
surements, finding a lack of sufficient material to account 
for most polluted stars. 

This paper returns to the question of the origin of the 
DZ spectral class. Historically, the DZ stars outnumbered 
their h ydrogen-rich counterparts stars by 25:1 |Dupuis et al.l 
1 1993b! ') until high resolution, spectroscopic searches with 
large telescopes uncovered a commensurate numbe r of DAZ 
stars (|Koester et al.ll2005al ; fZuckerman et al.ll2003T l. Today, 
the ratio of known DZ to DAZ stars is greater than 4:1, due 
almost exclusively to the Sloan Digital Sky Survey (SDSS; 
lEisenstein et al.ll2006T ). The ability to detect mild levels of 
heavy element pollution in these stars, combined with the 
relative lack of circums tellar material at both the DZ and 
cool DAZ white dwarfs l|Farihi et alj|2009h , makes the origin 
of the DZ stars challenging to constrain. 

It is argued here that the hundreds of newly identified, 
cool, helium-rich white dwarfs from the SDSS provide an 
excellent, spatially and kinematically unbiased sample upon 
which to test the ISM accretion hypothesis. These stars form 
a large statistical sample that favors an alternative to the 
accretion of Galactic gas and dust; as a whole the SDSS 
DZ stars do not show the spatial, kinematical, or elemen- 
tal abundance correlations expected from accretion episodes 
within the plane of the Galaxy. Furthermore, the SDSS DZ 
and DC (helium-rich, featureless) stars appear to be similar 
populations of polluted and non-polluted stars broadly shar- 
ing all relevant characteristics save photospheric calcium, 
and perhaps hydrogen. 



2 COOL HELIUM-RICH WHITE DWARFS 
FROM THE SDSS DR4 

The SDSS fourth data release (DR4; 

lAdelman-McCarthy et ail I2006T I contains a prodigious 
numb er of new white dwa rfs, in the neighborhood of 
6000 (jEisenstein et al.1 I2006T I. Excluding multiple entries 
for the same source, binaries or binary candidates, and 
uncertain classifications, there are 95 DZ and 142 DC 
stars wit hin the DR4 white dwarf catalog. Under further 
scrutiny, iDufour et al.l (|2007h has expanded the number of 
confirmed, unique DZ stars to 146, providing effective tem- 
peratures, calcium abundances, and photometric distances 
under the assumption of log [g (cms" 2 )] = 8.0. Of these 
288 cool, helium-rich white dwarfs in DR4, only one was 
known previously; Glll-54 or SPSS J080537.64+383212.4 
l|Dufour et al.ll2007l : [Gr"eensteinlll975l ). 

Perhaps surprisingly, the current number of DC and 
DZ stars in the SDSS are nearly identical. The DR4 cat- 
alog autofit temperatures for helium-rich stars are gener- 
ally unreliable near or below 10 000 K, as their helium at- 
mosphere models do n ot extend below this temperature 
|Eisenstein et al.1l2006T ). Owing to this fact, the DC star 
spectroscopic and photometric data were fitted using the 
same techniq u e and helium-rich (metal-free) models as in 
IDufour et all (|2007l ). yielding effective temperatures and 
photometric distances, again assuming log g = 8.0. All pho- 
tometric data were tak en from the DR4 white dwarf catalog 
l|Eisenstein et al.1 20061). with s tellar parameters for the DZ 
stars from IDufour et al.l (|2007T l and the DC model fits per- 
formed for this work. 

Spatial and kinematical data for these white dwarfs 
were obtained from t he SDSS seventh data release (DR7; 
lAbazaiian et al.ll2009T l. This latest data release includes the 
USNO -SDSS derived proper motio n catalog of Munn et al.l 



|2004l) and its recent amendments (|Munn et al.ll2008[) . with 
statistical errors around 3masyr _1 . For a handful of white 
dwarfs, proper motions were not available in the DR7 cata- 
log; for these stars measurements were taken from the S uper- 
COSMOS Scien ce Archive (SSA: iHamblv et al.ll200lf ), the 
LSPM catalog jLepine fc Sharal I2005T I . or were calculated 
for this work based on two or more positions on archival 
photographic plates and SDSS images separated by approx- 
imately 50 years. These proper motions are listed in Table 

m 



3 ABUNDANCE ANALYSIS 

The SDSS white dwarfs offer a deeper snapshot of the Solar 
neighborhood and beyond. Owing to the g > 22 AB mag- 
nitude limit of the survey, and its general avoidance of the 
Galactic plane, these white dwarfs typically represent dis- 
tant stars (tens to hundreds ofpc) at intermediate Galactic 
latitudes. Hence the SDSS DZ stars are an excellent probe 
of the local ISM (LI SM) both within and be yond the Local 
Bubble or Chimne y (IWelsh et al.ll 19991. 1 1994) . which extends 
roughly to 100 pc (|Redfield fc LinskvlBoOsf) . 

A number of relevant quantities were calculated from 
the DZ and DC white dwarf spatial and kinematical data: 
tangential speed utan, height above (or below) the Galactic 
mid-plane \z\, space velocities UVW (assuming zero radial 
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Table 1. White Dwarf Proper Motions Unavailable in SDSS DR7 



SDSS White Dwarf fi P.A. Sources 

(masyr - - 1 ) (deg) 

1,2,3 
1 
1 
3 
1 

1,2,3 
1 
1 
1 

2,3 
1 
3 
3 
1 



Sources: (1) SSA; (2) LSPM catalog; (3) this work. 

Table 2. Statistical Properties of SDSS DZ and DC White 
Dwarfs 



Parameter 


146 DZ Stars 


142 DC Stars 


T eff (K) 


8700 ± 1330 


9040 ± 1610 


9-2 


-0.29 ±0.23 


-0.18 ±29 


d (pc) 


196 ± 91 


206 ± 96 


\z\ (pc) 


140 ± 67 


148 ± 72 


ftan (kms -1 ) 


48 ±32 


52 ± 25 


U (kms" 1 ) 


2 ±35 


-4± 38 


V (kms' 1 ) 


-12 ± 28 


-12 ± 24 


W (kms" 1 ) 


1 ±25 


-1 ± 25 


\W\ (kms" 1 ) 


18 ± 17 


19 ± 15 


T (kms" 1 ) 


43 ± 29 


48 ±23 



Space velocities were calculated assuming zero radial velocity and 
corrected to the LSR (otherwise T = «tan). 

velocity), and total space velocity T 2 = U 2 + V 2 + W 2 . 
Space velocities were corrected to the Local Standard of Rest 
(LSR; iDehnen fc Binnevll 19981 ). with U positive toward the 
Galactic anticenter, V positive in the direction of Galactic 
rotation, and W positive toward the north Galactic pole. In 
the absence of correction to the LSR, the assumption of zero 
radial velocity is equivalent to T = utan- Tangential speeds 
were calculated from 

ftan (kms -1 ) = 4.7405 x d(pc) x ji (arcsec yr -1 ) (1) 

The Sun does not have zero velocity relative to the LSR; 
depending on viewing direction, the Solar motion con- 
tributes t; t an = 10 - 15 kms" 1 (|Dehnen fc Binnevl [l99l 
iMihalas fc Binnevlll98ll ). Table[2]lists the calculated param- 
eter means for the SDSS DZ and DC white dwarfs. 

3.1 A Glance at the Data 

As a first attempt to establish the relationship, if any, be- 
tween the DZ stars and the ISM, one can search for spa- 
tial or kinematical cor r elations with calcium abundance 
l|Zuckerman et al.l 120031 ; lAannestad et al.l Il993l ). Figure [1] 



plots the calcium abundances (throughout this paper [X/Y] 
= log [n(X)/n(Y)]), relative to helium, of the 146 DZ stars as 
a function of effective temperature, height above the Galac- 
tic mid-plane, and tangential speed. A low abundance cut- 
off at higher effective temperatures and larger distances is 
apparent in the upper and middle panels; these are observa- 
tional biases arising from higher atmospheric opacities and 
lower spectroscopic sensit ivity, respectively (jDufour et al.l 
120071 ; iKoester et aill2005al 'l . The classical conundrum of the 
DZ white dwarfs presents itself well in the upper panel; the 
majority of stars have no detected hydrogen, including those 
with the highest calcium abundances. 

The middle and lower panels of Figure [T] are notewor- 
thy because they fail to reveal a physical correlation be- 
tween accreted calcium and height above the Galactic disk 
or tangential speed. The Galactic gas and dust layer extends 
roughly 100 pc above the Galactic mid-plane, and yet some 
of most highly polluted DZ white dwarfs lie well above this 
region. The lower panel is relevant in the context of gravi- 
tationally driven accretion, expected if the DZ stars obtain 
their heavy elements while moving through the ISM. The 
lack of correlation between speed and the metal abundances 
argues, qualitatively but strongly, against fluid accretion of 
ISM. 

Figure [2] plots the masses of calcium contained in the 
convection zones of these metal-enriched stars using their 
abundances together with the conve ctive envelope masses for 
log g = 8 helium-rich white dwarfs (|Koesterj |2009). The fig- 
ure shows that the most highly polluted DZ stars currently 
harbor up to 10 22 g of calcium in their convection zones. This 
is a truly remarkable amount of a single heavy element, a 
mass in calcium alone that is roughly equivalent to the to- 
tal mass of a 200 km diameter asteroid(!). A more typical 
DZ star appears to contain around 10 20 g of calcium, still 
a prodigious amount, and correspo nds to a time-ave raged 
metal accretion rate of 2 x 10* g s _1 (|Farihi et al.ll2009l ), and 
a total accreted mass of metals of just under 10 22 g. 

Figure [3] plots the hydrogen abundances in 37 of 
the 146 DZ stars w here hydrogen is detected or inferred 
( Dufo ur et al. I l2007h . versus height above the Galactic mid- 
plane and tangential speed. No obvious pattern is seen, al- 
though there may be a higher density of DZA stars near the 
Galactic disk and perhaps also toward more modest speeds, 
but the former may be an observational bias due to dimin- 
ished spectroscopic sensitivity. In any case, one should ex- 
pect a correlation between these quantities if the hydrogen 
were accreted from ISM. 

As discussed by iDupuis et al.l lfl99 3a1. white dwarfs 
accreting Solar abundance matter within the ISM at 3 x 
lO^gs- 1 for 10 6 yr can accumulate the heavy element mass 
fractions observed in the convective envelopes of DZ stars, 
such as plotted in Figure [2] (however, this cannot account 
for their hydrogen abundances, discussed in detail below). 
Yet some care must be taken not to invoke such high ac- 
cretion rates without skepticism; while these rates are suf- 
ficient to account for the calcium data (that is the nature 
of the hypothesis), the physical plausibility is of equal im- 
portance. Below, the problem of interstellar accretion onto 
white dwarfs is reviewed in some detail and updated with 
old and new theoretical and empirical considerations. 
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Figure 1. Calcium abundances for the 146 DZ stars plotted ver- 
sus effective temperature, height above the Galactic mid-plane, 
and tangential speed. Those stars with detected hydrogen are 
plotted as filled circles, while those without are plotted with open 
circles. 



3.2 A Closer Look at the Interstellar Accretion 
Hypothesis: Hydrogen 

Both iKoesterl ( |l976l ) and IWesemaell (|l979l l were among 
the first to consider the issue of interstellar accretion onto 
helium-rich white dwarfs, both concluding that the atmo- 
spheres of these stars (both type DB and DC) were incom- 
patible with Bondi-Hoyle accretion of interstellar hydrogen. 
That is, in most cases a single interstellar cloud encounter 



Figure 2. Mass of calcium and hydrogen in the convective en- 
velopes of_the SDSS DZ stars; envelope masses were taken from 
IKoesterl <2009l) . A typical calcium mass in the convection zone of 
these stars is 6 X 10 19 g, roughly the same as that contained in 
the asteroid Minerva, whose mean diameter is 146 km and total 
mass is 3 X 10 21 g. 



lasting around 10 5 yr would suffic e to transform a DB (or 
DC) star into a DA (Kocstcr 1976). Therefore, the existence 
of more than 1000 hyd rogen-poor white dwarfs in the SDSS 
(jEisenstein et al.ll2006r ) argues strongly against the Bondi- 
Hoyle type accretion of interstellar hydrogen. 

Regardless of particulars relevant only to the metal- 
enriched white dwarf varieties, the lack of hydrogen in 
helium atmosphere white dwarfs is of great observational 
significance. High resolution and high signal-to-noise spec- 
troscopy reveals roughly equal numbers of DB stars where 
hydrogen is either undetected or detected but deficient by 
a few to several orders of magnitude relative to helium 
l|Voss et al.ll2007l ). The highest mass of hydrogen seen in a 
DB (metal-free) atmosphere is around 8x 10 22 g in a 12 000 K 
star, while the older and cooler counterparts studied here 
(the DZA stars) contain up to 4 x 10 g. This higher mass of 
atmospheric hydrogen would be accreted in u nder 10 6 yr at 
the int erstellar Bondi-Hoyle rate assumed by iDupuis et al.l 
(|l993ah . and hence this cannot be the correct picture due to 
the lack of such hydrogen masses in DB stars. 

Gravitational accretion of ISM particles by a star of 
mass M and radius R, in the supersonic regim e follows the 
mathematical form (jAlcock fc Illarionovlll98(l ) 
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(|Bondilll952h . yielding an accretion rate for interacting par- 
ticles 
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Figure 3. Hydrogen abundances for the 37 DZA stars plotted 
versus height above the Galactic mid-plane, and tangential speed. 



M = 7T (R A + R) Rsp 



(2) 



where Ra is the accretion radius, s — \/v 2 + c a 2 for relative 
stellar velocity v and ambient sound velocity c s , while poo 
is the unperturbed density of material being accreted. The 
accretion radius is defined as 



Ra 



2GM 



(3) 



and is often referred to the as the Bondi or Bondi-Hoyle 
radius. For all possible speeds considered here, Ra 3> R, 
and the accretion rate is given by 



M = 2-kGMRaRspo 



Me 



2-kGMRpo 



(4) 



(5) 



This is the Eddington rate ( Ed dingtonl 1 1926T ) . the accre- 
tion induced on non-interacting particles by the geometrical- 
gravitational cross section of the star as it travels through 
the ISM. 

Bondi-Hoyle theory (i.e. including gas pressure; Edgar 

demonstrates that the effective cross section in Equa- 
tion [4] ttRaR becomes tvRa 2 in the fluid dynamical limit 



This mass infall rate represents the maximal, idealized case 
where the mean free path of the particles is such that colli- 
sions are important and transverse momentum is effectively 
destroyed downstream from the star. It is also physically 
unrealistic in perhaps all situations excepting an ionized 
plasma, as it assumes no net angular momentum between 
the a ccreting star and its surrounding medium |Koesterl 
Il976h . and it c ertainly does not apply to neutral atoms or 
large particles (|Alcock fc Illarionovlll980h . The ratio of the 
Bondi-Hoyle to Eddington accretion rates is v 2 R/2GM for 
v ^> c, or arou nd 10 4 for typical white dwarf sizes and speeds 
l|Koesterlll976h . Table [3] lists typical densities and other rel- 
evant parameters for four fundamental types of ISM: molec- 
ular clouds, diffuse clouds, warm ionized, and hot ionized. 
Listed also are the expected high-end mass infall rates for 
white dwarfs moving through these regions following either 
Bondi-Hoyle (fluid) or Eddington (geometric) type accre- 
tion. 

The two-phase accr etion-diffusion scenario as laid out 
by (|Dupuis et al.lll99 3a), invokes 10 6 yr within a cloud (ac- 
cretion), and 5 x 10 7 yr between clouds (diffusion). For disk 
stars with Galactic orbits similar to the Sun, these timescales 
imply about five cloud encounters per 240 Myr orbit. The 
motions of interest are the relative motions between stars 
and the ISM, the latter which should be moving within 
the relatively low veloci ty spiral arms of the Galactic disk 
(|jahreifi fc Wielenlll997T l. 

If correct, this implies a typical helium-rich white dwarf 
moving at 50kms _1 relative to the LSR and (presumably) 
the ISM travels, on average, 240 pc within clouds during a 
single Galactic rotation. This timescale corresponds to the 
cooling age for a log g — 8.0, 14 500 K helium-rich white 
dwarf, and hence such a star should, according to this sce- 
nario, obtain up to 5 x 10 22 g of hydrogen via geometric ac- 
cretion, or [H/He] = —4.4. This value is commensurate with 
the highest hydrogen abundances observe d in DB stars, an d 
well above the lower limit of detectability l|Voss et al.ll2007T ). 

This comparison implies 1) the densities and corre- 
sponding Eddington accretion rates in Table are too high, 
or 2) cloud encounters last less than 10 6 yr or are less fre- 
quent than once per 5 x 10 7 yr, but is otherwise consistent 
with the geometric capture of hydrogen within the ISM, and 
inconsistent with fluid accretion. Ignoring the fact that all 
hydrogen will be ionized within the Bondi-Hoyle radius of 
the white dw arf, and hence should acc rete as a plasma at 
the fluid rate (lAlcock fc Illarionovll 19801) . it is clear that fluid 
accretion of hydrogen does not occur in helium atmosphere 
white dwarfs, if and when passing through dense ISM. 

3.3 Super-solar Calcium to Hydrogen Abundances 

The average DZ star in the sample has an effective tem- 
perature of 9000 K and a cooling age around 860 Myr 
|Fontaine et alj l200ll) . A hydrogen abundan ce of [H/He] 
= — 3.8 was detectable in all 146 DZ stars l|Dufour et al.l 
l2007h . equivalent to a bit less than 10 24 g of hydrogen at the 
average DZ temperature. In order to avoid such a hydrogen 



© 0000 RAS, MNRAS 000, 000-000 



6 J. Farihi et al. 



Table 3. Representative High ISM Accretion Rates for White Dwarfs 



ISM Type Gas Poo T c s M B H M Edd 

(cm" 3 ) (K) (kms" 1 ) (gs" 1 ) (gs" 1 ) 



Molecular Cloud 


H 2 


1000 


20 


0.3 


10 12 


10 8 


Diffuse Cloud 


H 


100 


80 


1 


10 11 


10 7 


Warm Ionized 


H+ 


1 


8000 


1 


10 9 


10 5 


Hot Ionized 


H+ 


0.01 


10 6 


100 


10 7 


10 3 



Accretion rates listed are order of magnitude estimates calculated for the 
listed parameters and v = 50 kms -1 , a speed typical of both the DZ and DC 
samples. 



mass within its outer envelope, a helium-rich white dwarf 
must eschew 1) fluid accretion in any region similar to the 
relatively low density LISM for 99.9% of its cooling age, and 
2) geometrical accretion in any molecular cloud, for around 
90% of its cooling age. Only 12 of the DZ stars or 8% of the 
sample have hydrogen abundances [H/He] —3.8, implying 
such avoidance has occurred for the remaining 134 DZ(A) 
stars. If the photospheric calcium originates in the ISM, then 
a process that accretes metal-rich dust grains while eschew- 
ing hydrogen (and other gaseous volatiles) is necessary to 
account for the DZ white dwarfs with little or no hydrogen. 

Historically, the super-solar [Ca/H] abundances of DZ 
white dwarfs were qualit atively accounted for by preferential 
accretion of dust grains (|Dupuis et al. I ll993bl ). A hypothesis 
to permit the incursion of dust grains yet repulse light gas is 
the propeller mechanism, w hereby charged ion species are re- 
pelled at the Alfven radius l| Wesemael fc Tru ran 1982). This 
model fails to explain a number of very hydrogen-deficient 
DZ-type stars, i ncluding GD 40 with [H/He] = —6.0, a white 
dwarf for which lFriedrich et al.l (|2004T ) failed to find evidence 
of magnetism sufficient to drive a propeller. Rather, in the 
same study the DZA star LHS 235, with 400 times more hy- 
drogen than GD 40 was found to have a 7 kG magnetic field 
- essentially the exact opposite of what might be expected 
if hydrogen was screened efficiently by a propeller. Perhaps 
the coup de grace for the propeller is the magnetic DZA star 
G165-7; with a 650 kG fiel d, [Ca/He] = -8.1, and [H/He] 
= -3.0 (|Dufour et al.ll2006l ). it has everything a 'fan' could 
hope for. But should this star have so much hydrogen? Ac- 
cording to theor y, its propeller failed arou nd 660 Myr ago 
as a 7500 K star (I Wesemael fc Truranll 19821 ). and could have 
amassed at most 2 x 10 24 g in 660 Myr of continuous accre- 
tion at the high Eddington rate, corresponding to [H/He] 
= —4.0. Clearly, this star was never a propeller. This mecha - 
nism is unsupported by observations |Friedrich et al.|[2004h . 
and is fund amentally at odd s with the existence DBAZ and 
DZA stars (|Voss et al.ll2007l ). 

The upper panel of Figure 2] plots the hydrogen ver- 
sus calcium abundances of all 37 DZA stars. The lack of 
any correlation between the two elements suggests either 
1) disparate origins, or 2) a common origin with intrinsic 
variations in [Ca/H]. Neither of these two possibilities is 
consistent with ISM accretion. The middle and lower pan- 
els of Figure 0] display the inferred [Ca/H] abundances and 
lower limits for the 37 DZA and 109 DZ stars. In both cases, 
the vast majority of these values are super-sol ar by at least 
an or der of magnitude (see also Figure 12 of iDufour et al.l 
l2007h . 



3.4 A Closer Look at the Interstellar Accretion 
Hypothesis: Metals 

Interstellar gas - both cold and warm - is metal-poor; heav- 
ier elements tend to b e locked up in grains, lea ving only light 
and volatile elemen ts dVenn fc Lambert! 1990). This fact was 
not considered by iDupuis et alj (|1993aT ). whose accretion- 
diffusion scenario explicitly assumes 1) accretion based on 
fluid rates for hydrogen also apply to heavy elements, and 
2) elements are accreted in Solar proportions. Above it was 
shown that fluid accretion rates for ISM accretion of hydro- 
gen over predict the observed abundances in DB stars by a 
few to several orders of magnitude. But more importantly 
for heavy element accretion is the fact that interstellar dust 
grains should not accrete at the fluid rate. 

Solid particles will follow trajectories independent of 
surrounding gas whether neutral, or charged and electri- 
cally or magnetically coupled to the gas (|Alcock fc Illarionovl 
1980). However, dust grains approaching the stellar surface 
will be separated into their constituent elements via sub- 
limation that proceeds at slightly different rates according 
to the particular chemical species, resulting in significant 
ingress prior to becoming gaseous. The gradual evapora- 
tion of approaching interstellar dust has two major conse- 
quences. First, gaseous heavy elements cannot accrete at the 
fluid rate, as the effective gravitational cross section is sig- 
nificantly reduced by the close approach necessary for their 
evaporation out of grains. Second, the process of sublima- 
tion will spatially and temporally separate heavy elements 
according to their specific evaporation temperatures, result- 
ing in preferentially higher accretion rates for volatile species 
such as ice mantles, as opposed to metallic atoms or refrac- 
tory material. This latter consequence is interesting because 
it predicts a relatively high volatile to refractory element 
accretion mixture; in stark contrast both to the observed 
abundances and to t hose expected from the accretion of cir- 
cumstellar material iSion et alj|l990ah . 

The radius at which dust grains are evaporated, 
i?ev, and become part of the ambient, interacting gas 
is substituted for R in Equation [2] which then becomes 
(|Alcock fc Illarionov|[l98cl ) 

Mgr = IT (R A + Rev) R CV Spoa (7) 



If the grains absorb and radiate energy as blackbodies, they 
will evaporate at temperature T ev at a distance from the 
star given by (|Chen fc Jura| [2001) 
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Figure 4. The upper panel plots the hydrogen versus calcium 
abundances for the 37 DZA stars, while the middle panel plots 
a histogram of their [Ca/H] ratios. The lower panel plots [Ca/H] 
ratio lower limits for the 109 DZ white dwarfs without detected 
hydrogen. The Solar abundance is [Ca/H] = —5.6. 



Taking a typical DZ effective temperature of 9000 K from 
Tableland T ev between 1500 and 2000 K, the resulting 
evaporation radius is between 10 and 18 white dwarf radii 
(or 0.13 to 0.23 R@). This is a few hundred times smaller 
than a typical Bondi-Hoyle radius, around 0.5 AU for a white 
dwarf, and hence heavy elements will accrete at a rate much 
closer to the Eddington rate. In this situation Ra S> Rev, 
and the resulting mass infall becomes 



All else being equal, this rate should be around 10 to 18 
times the Eddington rate for refractory-rich dust grains, 
such as inter stellar sil i cates. As pointed out more than thirty 
years ago bv lKoesterl l|l97rj ). the above analysis ignores any 
net angular momentum between the star and the surround- 
ing ISM, which will work to reduce the effectiveness of ac- 
cretion excepting where material is fully ionized, a situation 
unlikely to be realized for gaseous heavy elements. The phys- 
ical considerations discussed in this section have been, for 
the most part, overlooked in several major studies of DZ 
stars, and their absence has biased outcomes favoring the 
interstellar origin of met a ls in cool white dwarfs in general 
(IKoester fc Wilkenl l20od: IWolff et all |2002| ; iFriedrich et al.l 
120001 ; iDupuis et al.lll993bl ). 

Applying Equation[8]to the highest density, Solar abun- 
dance ISM (i.e. a dust to gas ratio of 1:100, [Ca/H] = —5.6, 
poo = 1000 gem 3 ), a white dwarf can accrete heavy elements 
from the ISM at up to 2 x 10 7 gs _1 . Over 10 6 yr of sus- 
tained accretion at this rate, a typical white dwarf would 
accumulate 7 x 10 18 g of calcium in its convective envelope; 
longer timescales are irrelevant because diffusion begins to 
compete with accretion. This quantity of calcium is able to 
account for a small minority of the observed envelope masses 
plotted in the upper panel of Figure [2] the average calcium 
mass is 10 times higher, with a few values 3 orders of magni- 
tude greater. This calculation assumes 1) the highest accre- 
tion rate for dust grains within dense ISM regions, and 2) 
stars that currently retain their peak calcium abundances to 
within an order of magnitude. Hence, the calculation should 
over predict the amount of calcium polluting the stellar con- 
vection zones by one to a few orders of magnitude. Therefore, 
if this analysis is valid to within an order of magnitude, it 
is all but certain that the photospheric calcium in DZ stars 
cannot have interstellar origins. 

The DAZ white dwarfs make this point even more strik- 
ing. For these polluted, hydrogen atmosphere stars with thin 
convection z ones and short metal diffu sion times relative to 
the DZ stars, iKoester k, Wilkenl (|200tJ ) infer continuous and 
ongoing mass accretion rates based on a highly probable, 
current steady state balance between accretion and diffu- 
sion. The inferred DAZ accretion rates for Solar abundance 
matter, up to 2 x 10 11 gs~ 1 , are qui te similar to the high 
rate assumed bv lDupuis et al.l (|l993al ) for the DZ stars. Yet 
no correlation is found between the implied DAZ accretion 
rates and the total space vel ocity of the stars as exp ected 
in the case of fluid accretion (|Koester fc Wilked[200tJ ). Fur- 
thermore, there are no known regions of enhanced LISM 
density in the vicinity of the DAZ stars dKilic fc Redfieldl 
l2007l ; IKoester fc WilkenlbOOfj ; IZuckerman et al.ll2003T l. and 
Table [3] reveals that accretion rates of 10 11 gs _1 are only 
even theoretically possible within relatively dense clouds, 
which are not seen. Lastly, there is firm evidence that the 
nearby DAZ G29-38 does not accrete matter at the fluid 
rate, but rather has an uppe r limit based on X -ray observa- 
tions that is 25 times lower l|jura et al,ll2009bh . 

There is no physical reason why helium and hydrogen 
atmosphere white dwarfs should become polluted via sep- 
arate mechanisms; in fact, Occam's razor argues against a 
bimodal theory. The fact that some examples of both va- 
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rieties of metal-contaminated stars are found with closel y 
orbiting dust is testament to this fact l|Farihi et al.l f2009). 
In this regard, the only real difference between the DZ and 
DAZ stars is that the latter group is currently accreting from 
its immediate surroundings, while the former group carries 
their metallic scars long after a pollution event has ended 
i|Koesterll2009l ). 



4 SPATIAL - KINEMATIC AL ANALYSIS 

4.1 Polluted White Dwarfs in the Local Bubble 

While iKoester fc W ilkcn (200fj) obtain capture rates up to 
10 11 gs" 1 in the warm, ionized LISM within the Local Bub- 
ble which assumes fluid accretion is valid for heavy elements, 
there are still serious problems with this picture when ap- 
plied to the DZ stars. The indi vidually identified parcel s 
of LISM span at most a few pc (|Redfield fc Lmsky||2008fh 
corresponding to several 10 4 yr of traversal by a typical DZ 
white dwarf. Such an encounter in a typical r egion of LISM 
with poo = 0.1 cm -3 (|Redfield fc Linskvll2000l ) leads to M B h 
rates of 10 s gs -1 , and would provide only 10 16 g of calcium 
pollution if accreting Solar abundance matter. While such 
modest an amount of metal can effectively pollute the rel- 
atively thin outer layers of DAZ stars, this paltry calcium 
mass requires 10 000 LISM encounters to account for the 
average calcium mass in the DZ star envelopes. Of course, 
this presumes fluid accretion is viable for heavy elements, 
which has been shown above to be physically implausible, 
thus making legion fluid infall episodes unlikely in the ex- 
treme. Therefore, it is safe to conclude that within the Local 
Bubble, ISM accretion is inconsistent with the phenomenon 
of DZ stars, which on average travel only 50 pc during a 
10 6 yr metal sinking timescale. 



4.2 ISM Accretion Above the Galactic Disk? 

Prior to the SDSS, there were few (if any) cool 
(T eS ^ 12 000 K) wh ite dwarfs known beyond 200 pc 
l|McCook fc Sionl Il999h - Generally a hindrance for both 
magnitude limited as well as proper motion surveys, such 
distant cool degenerates would have V > 18 mag, while 
100 km s _1 tangent ial speeds only y i eld fi < 100 mas yr -1 . 
The DZ sample of lAannestad et al.l fl993) was limited to 
within 50 pc of the Sun, but the SDSS DZ sample stretches 
out to 300 pc above the disk of the Milky Way (see Figure 
[TJ. While the single DZ outlier at \z\ > 400 pc is likely to 
be a higher than average mass white dwarf (i.e. its pho- 
tometric distance is overestimated), the sample as a whole 
contains 96 (66%) and 26 (18%) stars situated more than 
100 and 200 pc above the plane of the Galaxy. This fact 
alone perhaps provides an insurmountable challenge to the 
ISM accretion hypothesis. 

The sc ale height of molecular ga s as measured by CO is 
75 ± 25 pc (|Dame fc Thaddeuj[l985l 'l . while that of neutral, 
atomi c hydrogen is roughly 135 pc with substantial irregular- 
ities (|Lockman et al.ll 19861 ). Taking 100 pc as a fiducial scale 
height for interstellar gas and dust, and assuming the entire 
sample of DZ stars obtained their metals from this medium 
requires the stars: 1) are currently traveling through or away 



from the Galactic mid-plane, and 2) have not traveled sig- 
nificantly further than around 18 pc above the gas and dust 
layer, on average. This distance constraint corresponds to a 
star traveling at the sample mean \W\ = 18kms _1 speed 
(see Table [2]) over a typical Myr diffusion timescale. 

However, these requirements do not match the obser- 
vations. First, more than half the sample (81 stars or 58%) 
currently sit more than 120 pc above the disk of the Galaxy. 
Second, and rather damningly, nearly half of the sample is 
moving back into the disk rather than away. The upper panel 
of Figure [5] reveals 67 polluted white dwarfs that are cur- 
rently moving back toward the Galactic mid-plane. Many of 
these stars are among the most highly polluted in the sam- 
ple, and many have \z\ > 100 pc. The middle panel of the 
same figure shows the height above the disk as a function 
of total space velocity; surprisingly some of the fastest stars 
are also high above the plane. These results strongly conflict 
with expectations from ISM accretion models. 

In order to estimate how much time high |z| DZ stars 
spend above the spiral arms and interwoven disk material, 
one can use the Solar motion as a first estimate. Depend- 
ing on the Galactic model potential, and whether it in- 
cludes dark disk m atter, the Sun oscil l ates about z = 
every 60 to 90 Myr jMatese et all 1 19951 ; IBahcall fc Bahcalll 
1 19851 ; iMihalas fc Binnevlll98ll '), the former and latter val- 
ues derived using models containing substant i al and no dark 
disk matter, respectively. iKuiiken fc Gilmorel (| 19891 ) provide 
VK m ax and z ma x data for 100 pc ^ z ^ 2000 pc based on a 
sample of K dwarfs toward the South Galactic Pole and 
a Galactic model which contains no dark disk component. 
Using a simple harmonic oscillator model, a period can be 
computed from these maximum heights and speeds for stars 
within several hundred pc of the disk. The period of a simple 
oscillator is 

P z =2-kJi— (10) 

V m 

where 

Kz = m (^L\ 2 (ii) 

V Zmax / 

In reality, the period is a function of z due to the disk 
component of the Galaxy, which acts as a linear term 
in the vertical gravitational potential, while the spherical 
component (including the inferred dark halo) acts as a 
qua dratic term, equival ent by itself to a harmonic poten- 
tial (IBahcall et al.ll 19991) ■ U s ing th e above equations and the 
data of IKuiiken fc Gilmorel 0989) for stars within 500 pc of 
the plane, the mean period is P z — 80±10 Myr and compares 
favorably to rigorously derived values for the Sun. Taking 
P z = 80 Myr, one can then use the same simple model to 
calculate the length of time a DZ star with current velocity 
W and height z has spent outside the 100 pc Galactic gas 
and dust layer. 

The lower panel of Figure [5] plots a histogram of the 
outward- and inward- bound DZ white dwarfs which are 
currently at least 100 pc above the Galactic mid-plane. All 
96 of these stars have spent more than 3 Myr in a region 
of space where the accretion of interstellar matter can be 
likely ruled out, 91 of these metal-polluted stars have been 
out of the Galactic disk for over 6 Myr. As can be seen by 
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If the majority of these DZ white dwarfs have traveled 
t ° away from where they were polluted for a few to several 

. o ° diffusion timescales, the extant photospheric metals repre- 
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Figure 5. The upper panel reveals 67 DZ stars moving back to- 
ward the Galactic mid-plane, many among the most highly metal- 
polluted stars in the sample. The middle panel shows the total 
space velocity as a function of height above the Galactic disk, 
with many stars exhibiting high values along both axes. The lower 
panel is a histogram of Myr spent at |z| > 100 pc for the 96 DZ 
stars currently located in this region. 



reviewing Figure [T] these high \z\ DZ stars are not anoma- 
lous in terms of their abundances, but rather mundane and 
truly representative of the sample at large. If anything these 
white dwarfs have higher abundances than the lower \z\ stars 
(but this is likely a bias resulting from spectroscopic calcium 
detection sensitivity) . 



4.3 ISM Accretion At High Speed? 

Throughout the preceding sections, a nominal 50kms _1 
tangential speed was used for relevant calculations; this 
value being somewhat friendly to ISM accretion models and 
also the mean among the 146 DZ stars in the sample. How- 
ever, there are 13 stars with Wtan ^ 100 km s -1 (and total 
space velocities, T > 100 km s" 1 since their radial veloc- 
ities are unknown and assumed to be zero), which would 
imply a decrease in the Table fluid accretion rates by an 
order of magnitude. There appears to be no correlation be- 
tween these speeds and the accreted calcium abundances in 
the DZ stars (see Figure [T]), and hence there seems little 
or no chance that the s e sta rs have accreted in this man- 
ner. lKoester& Wilkenl (2006) finds a similar lack of correla- 
tion among the known DAZ stars; a more striking statement 
given their necessary ongoing accretion. 

While the dependence of geometrical accretion on ve- 
locity is milder by comparison, for ISM accretion scenarios 
it is reasonable to expect a mild correlation should exist be- 
tween accreted mass and velocity, yet no such correlation is 
seen in Figures [T] or [3J In the middle panel of Figure [5] one 
would expect to find milder velocities necessary to capture 
more ISM material as stars move away from the Galactic gas 
and dust layer and into regions where the ISM is tenuous at 
best, yet this pattern is not observed. 



5 THE SDSS DZ AND DC WHITE DWARFS 
COMPARED 

If the ISM accretion hypothesis is correct, then the only dif- 
ference between a DZ and a DC star is its recent (within 
a few Myr) spatial and kinematical history. That is, dur- 
ing that time a DZ star has recently exited a relatively 
dense patch of ISM, while a DC star has not (the con- 
tinuous, warm-phase ISM accretion scenario forwarded by 
iKoester fc Wilkenl |2006j ) does not work for the DZ stars, as 
demonstrated above). Is this the correct picture? 

There are a similar number of DZ and DC white dwarfs 
in the SDSS and a comparison of their model effective tem- 
peratures, g — z colors (the longest baseline possible which 
is unaffected by calcium H and K absorption), Galactic 
heights, and tangential speeds is shown in Figure [6] There 
appear to be a few modest differences between the sample 
distributions; a small group of relatively lower temperatures, 
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Figure 6. Effective temperature, g — z color, Galactic height, and tangential speed histograms for the DZ and DC white dwarfs. For the 
most part, these two spectral classes appear to represent the same population of helium-rich white dwarfs. 



and a flatter color distribution are found among the DC 
stars, while the zj and n tan distributions are fairly similar. 
Table [2] demonstrates the two spectral classes have mean 
parameters which agree rather well despite the relatively 
large standard deviations; their average effective tempera- 
tures, distances, Galactic heights, and tangential speeds all 
agree to within 4% to 8%. 

Figure [7] plots the UVW space velocities of the DZ 
and DC samples (assuming u ra d = 0), supporting the idea 
they are both basically thin disk populations, each with sev- 
eral possible thick disk members. Table [2] shows the two 
groups share the same mean space velocities and velocity 
dispersions, the l atter being the most important indicator of 
kinem atical age l|jahreifi fe W iden 1997; Mihalas " fc Binnevl 
Il98ll ). Together these analyses suggest DZ and DC white 
dwarfs belong to the same population of helium atmosphere, 
thin (and a few thick) disk stars. 



5.1 Convergent and Divergent Pairs of DZ and 
DC White Dwarfs 

Figure [S] projects the Galactic longitude and latitude of all 
the DZ and DC stars onto the plane of the sky. While it 
is true that the SDSS observes only certain areas of the 
sky, and hence one can a priori expect these two groups to 
occupy overlapping regions, the degree of overlap between 



the DZ and DC white dwarfs is total; there does not exist 
any region where DZ stars are found in higher concentrations 
than their DC counterparts. However, this type of projection 
is very dense on paper while the actual sky is rather vast, 
and it necessarily ignores the third dimension, namely the 
distance from the Sun. Therefore a pertinent question to 
ask is whether any of these stars are pc-scale neighbors. A 
search was conducted around each DZ and DC star for other 
sample stars within a r = 10 pc radius, and the results are 
listed in Table H 



All together 18 pairs of stars were identified in this 
manner, three of which lie together within a 5pc region 
of space. There are nine DZ-DC, 4 DZ-DZ, and 5 DC-DC 
pairs; exactly half the pairs are of mixed spectral type. Of 
the DZ-DC pairs, there are three where the total space ve- 
locities agree to within 12kms _1 , pairs #1, 6, and 10. In the 
latter case, each space velocity component by itself agrees 
within 3kms (!). These three pairs of stars have traveled 
the last 10 6 yr in adjacent regions of interstellar space, yet 
each component star has evolved distinctly from its partner 
- one metal-polluted and the other not. Furthermore, there 
are two DZ-DZ pairs whose total space velocities differ by 
more than lOOkrns -1 , thus implying their current spatial 
convergence is in contrast to their divergent interstellar his- 
tory over the last 10 6 yr. Therefore, these two DZ-DZ pairs 
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Figure 7. UVW space velocities for the DZ and DC white 
dwarfs calculated assuming f ra d=0. The ellipses shown in the 
plots are the 1 and 2a contours for old, metal-poor disk stars 
from lBeers et al.l |200Cft . 

obtained their metals while traversing distinct rather than 
adjacent regions of the Galaxy. 

Both the DZ-DC pairs with convergent spatial histories, 
as well as the DZ-DZ pairs with divergent orbital pasts are 
strong counterexamples to expectations for accretion within 
dense interstellar clouds of sizes several pc or larger. This 
is not surprising in the sense that such clouds do not exist 
within the Local Bubble, but again necessitates the exis- 
tence of plentiful, smaller ISM patches to account for 1) 
many of the pairs in Tab le [H and 2) the DAZ spectral class 
l|Koester fc Wil kcn 2006). As discussed above, the scenario 
of small (undetectable) cloud encounters may work for the 
DAZ stars, but not for the DZ class; it requires thousands 
of such encounters to accumulate the large mass of metals 
currently residing in their convection zones. 

5.2 The Frequency of Metal-Free Counterparts to 
DZA Stars 

A white dwarf accreting Solar ratio abundances at the high 
Eddington rate for 10 6 yr will obtain about 3 x 10 21 g of hy- 
drogen, and it is tempting to imagine the hydrogen masses 
in Figure |2] repres ent an accumulation of such encounters. 
lJura et aU (|2009bl ) plots a similar figure but in cluding hy- 
drogen masses in warmer DBA white dwarfs from lVoss et al.l 




o 
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Figure 8. Projection of the DZ and DC Galactic coordinates. 



(2007). Generally speaking, the hydrogen envelope masses in 
DBA stars are one to a few orders of magnitude smaller than 
those in the DZA, making it plausible that as helium-rich 
stars cool and accrete ISM gas, the hydrogen accumulates. 

If helium-rich stars gradually accumulate hydrogen 
from the ISM as they cool, there should exist a signifi- 
cant population of cool, helium-rich white dwarfs that reveal 
themselves as DA stars, which has not yet been observed 
|Zuckerman et al.ll20o3 : Bergeron et al.ll200ll ). Accreted hy- 
drogen should perpetually float in a helium-rich star, while 
any metals eventually sink below detectability; this expec- 
tation is corroborated in high-resolution spectroscopy of DB 
stars, where the metal -free DBA outnumber the DBAZ stars 
15:1 l|Voss et al.l 120071 ). While this ratio should decrease in 
the cooler DZA temperature range as diffusion timescales 
and atmospheric opacity favor a more frequent detection of 
calcium, the relative number helium-rich DA and DZA stars 
should remain larger than unity for all scenarios where met- 
als and hydrogen have distinct origins. This arises because 
separate sources for hydrogen and metals, in helium atmo- 
sphere white dwarfs, imply independent probabilities, and 
being polluted by both will be the least probable. Therefore, 
if one assumes that metals result from accreted asteroids 
but under-abundant hydrogen in a cool, helium atmosphere 
is either 1) interstellar, 2) primordial, or 3) the result of 
dredge up and subsequent mixing, a larger than unity ratio 
of metal-free counterparts to the DZA stars is certain. 

However, if the photospheric metals and hydrogen both 
have a circumstellar origin, and are semi-continuously de- 
livered to the photosphere on relevant timescales, the ra- 
tio of helium-rich DA to DZA stars can potentially be less 
than unity. There are three mixed helium-hydrogen atmo- 
sphere stars among the 1 52 cool white dwarfs t horoughly 
studied and modeled by iBergeron et al.l |2001); two are 
DZA stars('), while the third has atmospheric carbon (where 
dredge up of interior carbon and helium in a previously 
hydrogen-rich atmosphere is a possibility). While this ra- 
tio of 1:2 (or 0:2) results from small number statistics, it 
may nonetheless be telling. The significant populations of 
DC and DZ stars, coupled with the relative lack of cool, 
helium-rich DA white dwarfs supports a circumstellar ori- 
gin for hydrogen in DZA stars. 

Furthermore, Occam's razor favors a single mechanism 
that can account for all the data simultaneously; hydro- 
gen delivered together with heavy ele ments in w a ter-rich 
minor planets is one such possibility. lJura et all (|2009bh 
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Table 4. Pairs within 10 pc of one another 



Pair 


Stari 


Star2 


SpTi 


SpT 2 


Ar 


AU 


AV 


AW 












(PC) 


(kms- 1 ) 


(kms" 1 ) 


(kms -1 ) 


1 


J020132. 24-003932.0 


J015604.19-010029.3 


DZ 


DC 


3.4 


6.1 


4.3 


1.8 


2 


J084828. 00+521422.5 


J083736.58+542758.3 


DZ 


DC 


5.4 


49.3 


57.6 


42.9 


3 


J084828. 00+521422.5 


J085133.95+542601.2 


DZ 


DC 


5.2 


20.5 


14.2 


22.1 


4 


J084828. 00+521422.5 


J090729.10+513805.6 


DZ 


DC 


9.2 


19.6 


7.1 


20.6 


5 


J093423. 17+082225.3 


J093819.84+071151.6 


DZ 


DC 


3.1 


7.8 


21.8 


13.1 


6 


J104915. 06-000706.2 


J105632. 20-001041. 3 


DZ 


DC 


8.9 


9.3 


7.3 


2.4 


7 


J123455.96-033047.1 


J124006.36-003700.9 


DZ 


DZ 


7.6 


2.1 


21.3 


12.7 


8 


J124006. 36-003700.9 


J124333.86+031737.1 


DZ 


DC 


8.4 


35.3 


0.6 


7.4 


g 


T1 31 33fi 964-573800 5 


11 31 750 20-1-600532 9 


DZ 


DC 


6.7 


9.1 


115.1 


55.6 


10 


J135137.07+613607.0 


J135137.07+613607.0 


DZ 


DC 


6.0 


2.4 


1.5 


1.6 


11 


J144022. 52-023222. 2 
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offer this hypothesis to explain, with a single pollution 
event, the collective facts available on the spectacularly 
metal-enriched, mixed helium-hydrogen atmosphere white 
dwarf GD 362. This disk-polluted, helium-rich star has an 
anomalously high hydr ogen abundance at [H/He] = —1.1 
l|Zuckerman et al.ll2007t). so much that it was originally clas- 
sified as DAZ l|Gianninas et al.1 120041 ) . GD 16 is a similar 
case w ith significant atmospheric hydrogen |Koester et al.l 
2005b) and a dust disk po lluting its helium-dominated atmo- 
sphere (jFarihi et al.ll20oljl ). while the two remaining helium- 
rich white dwarfs with disks (GD 40 and Ton 345) have little 
or no hydrogen |jura et al.ll2009bl ). Therefore, the pattern 
of hydrogen abundances in DZ stars is likely a reflection of 
the diversity of water content in extrasolar planetesimals. 



6 OUTSTANDING ISSUES 

Owing to the fact that radial velocities are not reliably mea- 
surable from the calcium lines in the SDSS spectra of the DZ 
stars, some of the derived kinematical quantities are mod- 
estly uncertain. Fortunately, the galactic latitudes of the 146 
stars are bound within 25° < |6| < 75°, and potential biases 
are mild at worst. For example, although the calculation of 
W velocity here is missing the third dimensional ingredient, 
which could potentially alter its sign (and hence direction), 
a non-zero radial velocity will tend to increase their speeds 
in any given direction [T 2 = U 2 + V 2 + W 2 = v tan 2 +u rad 2 ]. 
Moreover, the focus of the study has been statistical in na- 
ture an d the global properties of the sample a re unlikely to 
change l|Pauli et al.ll2003l ; ISilvestri et ai]|2002l ) 

GAIA will eventually obtain distances to, and excel- 
lent proper motions for, most if not all the DZ and DC 
white dwarfs, the latter stars being immune to radial velocity 
measurements by nature (until more powerful instruments 
resolve any potential lines in their spectra). Empirical dis- 
tance determinations will remove any bias introduced with 
the assumption of log g = 8.0, but again this is unlikely to 



be pr oblematic from a statistical point of view (Kepi er et al.1 
l2007t) . 

Perhaps more important is the current understanding 
of how interstellar accretion occurs; this ignorance alone 
can sweep up many uncertainties, just as passing stars may 
sweep up ISM, but this situation is not ideal, and empirical 
constraints are superior. Although the propeller mechanism 
hypoth esis may be currently broken, or minimally in need of 
repair (|Friedrich et al.|[2004l ). this does not rule out the log- 
ical possibility of preferential accretion for heavy elements 
(grains) versus gas while a star traverses a patch of ISM. 
However, there are two worthwhile considerations for any 
such hypothesis. 

First, based on reasonable physics in the absence of 
propeller-like deterrents, the infall of gas onto a star via the 
ISM should proceed at a rate much higher than the analo- 
gous accretion of grains (§3.4). Two additional factors should 
further increase the accreted volatile-to-refractory element 
ratio: 1) the intrinsic gas-to-dust ratio in the ISM is around 
100:1, and 2) interstellar grain species should be separated 
further in the accre tion process owing to vary ing evapora- 
tion temperatures (|Alcock fc Illarionovl Il980l ). At present 
there does not exist a large number of cool white dwarf 
spectra exhibiting sufficient elements to confidently evalu- 
ate the volatile-to-refractory element ratios. Obtaining such 
a database would be highly valuable. 

Second, there is a class of main-sequence star which is 
thought to accrete ISM and become polluted in the oppo- 
site sense to contaminated white dwarfs. The A Bootis stars 
are a class of late B to early F stars that show nearly So- 
lar abundances of light elemen ts but are deficient in met- 
als jBaschek fc Slettebak|[T98i ). It has been suggested that 
their abundanc e ratios mimic the meta l-depleted pattern of 
interstellar gas (|Venn fc Lambertlll990h , and that accretion- 
diffusion of t his volatile-enriched gas i s consistent with th e 
observations (|Kamp fc Paunzenl [20021 : ICharbonneaulll99ll ). 
Since A Bootis stars are generally A-type dwarfs where con- 
vection has become important, they should share some ba- 
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sic characteristics with cool white dwarfs in the appropri- 
ate effective temperature range. If the A Bootis stars are 
polluted by interstellar, metal-poor gas, then any ISM ac- 
cretion theory for white dwarfs should be able to account 
for both classes of stars, and explain specifically why pref- 
erential accretion would favor heavy elements in the case of 
white dwarfs, contrary to reasonable physical expectations. 



7 CONCLUSIONS 

The SDSS currently contains a large number of cool, metal- 
polluted, helium atmosphere white dwarfs over a relatively 
large volume of the local Galaxy. These DZ stars appear to 
be, on average, atmospherically, spatially, and kinematically 
similar to a comparable number of DC white dwarfs from 
the same catalog, at least superficially. A search for corre- 
lations between calcium and hydrogen abundances, spatial 
and kinematical properties among the DZ stars yields neg- 
ative results, excepting a tendency for increased hydrogen 
mass with cooling age in DZA stars. Specifically, the metal 
abundances of 146 DZ stars shows no connection with speed 
relative to the LSR, nor with height above the Galactic disk. 

The super-solar [Ca/H] ratios and lower limits are as ex- 
pected for the DZ spectral class, but nonetheless difficult to 
understand as accretion of matter which is over 97% hydro- 
gen and helium. The substantial fraction of the DZ sample 
which sits at \z\ > 100 pc is difficult to reconcile with accre- 
tion of ISM, perhaps impossible. The total space velocities 
of these high \z\ stars are no slower than their counterparts 
within the Galactic gas and dust layer, and hence no more 
amenable to accretion of rogue ISM patches, which are un- 
likely to exist in any case. These 96 stars have spent 3 to 
40 Myr above the 100 pc half- width scale height of ISM, with 
91 spending at least 6 Myr, the latter time period sufficient 
to drive their metal abundances down by a factor 400. 

The calcium masses currently residing in the convective 
regions of the DZ stars are near 10 20 g on average. If one as- 
sumes the calcium rep resents 1.6% of the total mass (equiva- 
lent to the bulk Earth; lAllegre et al.lll995h of heavy elements 
in the stellar envelope, then a typical DZ star contains a 
mass equivalent to a 150 km asteroid. This is the minimum 
mass of accreted heavy elements. If a large percentage of 
DZ stars have retained metals from a past encounter within 
a few to several Myr diffusion timescales, then the present 
envelope masses of heavy elements represent only a fraction 
of the accreted masses and invoking total captured masses 
as large as the Moon. 

Simply taking the 146 DZ stars out of 4193 DR4 white 
dwarfs with T e s < 12 000 K, one can say that at least 3.5% 
of white dwarfs appear to be polluted by circumstellar mat- 
ter, the remains of rocky planetary systems. This trans- 
lates directly into a similar lower limit for the formation 
of terrestrial planets at the main-sequence progenitors of 
white dwarfs, primarily A- and F-type stars of intermediate 
mass. While this fractio n is likely to be significantly higher 
(IZuckerman et al. 1 120031 ). it is difficult to quantify without 
a commensurate examination of all the cool SDSS white 
dwarfs, which is beyond the scope of this paper. The ap- 
pearance of hydrogen in DZA stars suggests a common origin 
for both heavy elements and hydrogen, and indicates DZA 
star pollution by water-rich minor planets may be semi- 



continuous on Myr timescales. This picture can be tested 
by searching for sensitive features of oxygen absorption in 
hydrogen- and metal-polluted white dwarfs. 
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